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Abstract 

The aim of this work is to study the spectral statistics of the asymmetric 
rotor model (triaxial rigid rotator). The asymmetric top is classically inte- 
grable and, according to the Berry- Tabor theory, its spectral statistics should 
be Poissonian. Surprisingly, our numerical results show that the nearest neigh- 
bor spacing distribution P(s) and the spectral rigidity Aa(L) do not follow 
Poisson statistics. In particular, P(s) shows a sharp peak at s = 1 while 
Az(L) for small values of L follows the Poissonian predictions and asymptot- 
ically it shows large fluctuations around its mean value. Finally, we analyze 
the information entropy, which shows a dissolution of quantum numbers by 
breaking the axial symmetry of the rigid rotator. 
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I. INTRODUCTION 



In the semiclassical limit [1,2] there is a clear connection between the behavior of classical 
systems (regular or chaotic) and the corresponding quantal ones. For quantal systems, 
corresponding to classical regular systems, the spectral statistics (-P(s) and A 3 (L)) follow 
the Poisson ensemble, while for systems corresponding to chaotic ones the Wigner ensemble 
is followed (see, for example, [3] and references therein). 

Nevertheless, some exceptions are known. The most famous case is perhaps the harmonic 
oscillator one, discussed in great detail in references [4,5]. It has been also found that low- 
energy spectral statistics of higher-dimensional separable Hamiltonian systems can show the 
level repulsion typical for chaotic systems. Especially critical in this sense are systems close 
to the harmonic oscillators and to rectangular wells [6]. 

The aim of this paper is to discuss another pathological case: the classically integrable 
triaxial rotator model (see, for instance, ref. [7]). Incidentally, this model has been used 
very often in the description of the low-lying states of the even-even atomic nuclei [8]. 

The asymmetric top described by the rotor model is a classically integrable system, but 
an analytical formula, as a function of quantum numbers, for its energy spectrum is not 
known. Nevertheless, numerical results can be obtained. By following the Landau approach 
[9], the Hamiltonian operator is split into 4 submatrices, corresponding to different symmetry 
classes. Each truncated submatrix is numerically diagonalized. Finally, the nearest-neighbor 
spacing distribution P(s) and the spectral rigidity A 3 (L) are calculated. Surprisingly, the 
spectral statistics of energy levels do not follow the predictions of Poisson statistics. 



Let us consider a system of coordinates with axes along the three principal axes of intertia 
of the top, and rotating with it. The classical Hamiltonian H of the top is given by 
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where J = (Ji, J2, J3) is the angular momentum of the rotation and a — 1/h, b — l/h, 
c = I/J3 are three parameters such that I±, I 2 and J 3 are the principal momenta of intertia 
of the top. The Hamiltonian is classically integrable and its action variables are precisely 
the three components J s , s = 1, 2, 3, of the angular momentum [7]. 

The quantum Hamiltonian H is obtained by replacing the components of the angular 
momentum, in the classical expression of the energy, by the corresponding quantum oper- 
ators Ji, J2 and J3. The commutation rules for the operators of the angular momentum 
components in the rotating system of coordinates are given by 

J r J s Js'Jr — ^rst Jt 1 (2) 

where e rst is the Ricci tensor and r,s,t = 1,2,3. Note that these commutation rules differ 
from those in the fixed system in the sign on the right-hand side [9]. 

As usual, the two operators J 2 = Jf + Jf + Jf and J 3 are simultaneously diagonalized on 
the basis of eigenstates | J, k) with integer eigenvalues J and k (k = —J,—J + 1,...,J — 1,J), 
respectively. The non-zero matrix elements of the quantum Hamiltonian H in the basis 
\J,k) are given by 

(J, k\H\J, k) = ^ (a + b){J{J + 1) - k 2 ) + — ck 2 , (3) 

(J, k\H\J, k + 2) = ( J, k + 2\H\J, k) = 

= —(a - b) y/(J - k)(J - k - 1)( J + k + 1)( J + k + 2) . (4) 

The quantum Hamiltonian H has matrix elements only for transitions with k — > k or k ± 2. 
The absence of matrix elements for transitions between states with even and odd k has the 
result that the matrix of degree 2 J + 1 is the direct product of two matrices of degrees J 
and J + 1. One of these contains matrix elements for transitions between states with even 
k, and the other contains those for transitions between states with odd k [9]. 
It is useful to introduce a new basis, given by 



\J,k;S) = -j={\J,k) + \J,-k)) , \J,0,S) = \J,0) , 



\J,k;A) = -j=(\J,k)-\J,-k)) , k^O. (5) 

By using this new basis, the total Hamiltonian matrix is decomposed in the direct product 
of 4 submatrices by considering the parity of the quantum number k: even (E) or odd (O), 
and the symmetry of the state: symmetric (S) or anti-symmetric (A). So the submatrices 
are labelled as follow: (E,S), (E,A), (0,S), (0,A). These are the classes of symmetry of the 
system. In Table 1 we show the dimension of each submatrix for a fixed J. 

The matrix elements of the Hamiltonian H in the new basis, with respect to the old 



basis, are given by 

(J,k,S\H\J,k,S) = (J,k,A\H\J,k,A) = (J,k\H\J,k) , k^l (6) 

(J, 1, S\H\J, 1, S) = (J, 1\H\J, 1) + (J, 1\H\J, -1) (7) 

(J, 1,A\H\J, 1,A) — (J, 1\H\J, 1) - (J, 1\H\J, -1) (8) 

(J,k,S\H\J,k + 2,S) = (J,k,A\H\J,k + 2,A) = (J,k\H\J,k + 2) , k ^ (9) 

(J,0,S\H\J,2,S) = V2(J,0\H\J,2) , fc ^ (10) 



We calculate the eigenvalues of each submatrix for different values of J using a fast imple- 
mentation, in double precision, of the Lanczos algorithm with a LAPAC code [10]. In Figure 
1 we plot the density of levels p(E) of each submatrix of H and J = 1000. The results show 
that the density of levels is practically the same for the four classes. p(E) displays a high 
peak at the left-center of the energy interval and a long tail for large energy values. 

III. SPECTRAL STATISTICS 

As previously discussed, according to the Berry- Tabor theory [11,12], given a classical 
integrable Hamiltonian that, written in action variables J r , satisfies the condition 
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(11) 



8J r dJ s 



then, in the semiclassical limit, its spectral statistics should follow the Poisson statistics. 
Note that a system of linear harmonic oscillators, whose Hamiltonian is given by H = uj • I, 
does not satisfy the previous condition. In fact, a system of linear harmonic oscillators is 
integrable but it does not follow Poissonian statistics [4,5]. 

The triaxial rigid rotator is integrable and satisfies the Berry- Tabor condition (11). Thus, 
one expects that the spectral statistics of the quantized rigid rotator should be Poissonian. 
We shall show that is not the case. 

In general, various statistics may be used to show the local correlations of the energy 
levels but the most used spectral statistics are Pis) and A 3 (L). P(s) is the distribution 
of nearest-neighbor spacings Sj = (E i+ i — Ei) of the unfolded levels Ei. It is obtained by 
accumulating the number of spacings that lie within the bin (s, s + As) and then normalizing 
P(s) to unit. As shown by Berry and Tabor [11,12], for quantum systems whose classical 
analogs are integrable, P(s) is expected to follow the Poisson distribution 



The statistic A 3 (L) is defined, for a fixed interval (—L/2, L/2), as the least-square deviation 
of the staircase function N(E) from the best straight line fitting it: 



where N(E) is the number of levels between E and zero for positive energy, between — E and 
zero for negative energy. The A 3 (L) statistic provides a measure of the degree of rigidity of 
the spectrum: for a given interval L, the smaller A 3 (L) is, the stronger is the rigidity, signi- 
fying the long-range correlations between levels. For this statistic the Poissonian prediction 
is 



P(s) = exp (— s) . 



(12) 
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It is useful to remember that Berry, on the basis of the Gutwiller semiclassical formula for the 
density of states, has shown that A 3 (L) deviates from the universal Poissonian predictions 
for large L: A 3 (L) should saturate to an asymptotic value performing damped oscillations 
[13]. 

In Figure 2 the spectral statistic P(s) is plotted for the four submatrices of H and 
J = 1000. Note that the level spectrum is mapped into unfolded levels with quasi-uniform 
level density by using a standard procedure described in [14] . As expected from the previous 
analysis of density of levels, P(s) is practically the same for the four classes of symmetry. 
Moreover, P(s) has a pathological behavior: a peak near s = 1 and nothing elsewhere. 
Compared to P(s), the spectral rigidity A 3 (L) is less pathological. As shown in Figure 3, 
A 3 (L) follows quite well the Poisson prediction A 3 (L) = L/15 for small L but for larger 
values of L it gets a constant mean value with fluctuations around this mean value. These 
fluctuations becomes very large by increasing L, in contrast with the Berry prediction [13]. 

The behavior of the density of levels p(E) and of the spectral statistics P(s) and A 3 (L) 
does not change by changing the matrix dimension, namely the quantum number J. In 
Figure 4 we plot the density of levels and the spectral statistics for J = 2000 and J = 4000. 

The results shown in the first four figures have been obtained with a = 1, b = \/2 and 
\/5, in such a way that the rotor is triaxial. It is interesting to see what happens if one 
changes the deformation parameters a, b and c, studying the transition from axial symmetry 
to triaxial symmetry. To do so, we take fixed b and c and modify a. In Figure 5 we plot the 
density of levels p(E) for six values of a ranging from a = a/2 to a = 0. The case a = b = y/2 
correspond to the axial symmetric one. The density of levels p(E) is strongly modified by 
changing the parameter a, i.e. breaking the axial symmetry, but the spectral statistics are 
not, as shown by Figure 6 for the P(s) distribution. 

To conclude this section, we discuss another statistical quantity that has been proposed 
to study quantum chaos: the information entropy S(E) of the eigenvector \E > associated 
to the eigenvalue E of the the Hamiltonian operator H [15]. Given a generic basis set {\i >}, 
the eigenvector \E > can be written as: 
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\E>= J £c i \i>, (14) 

i 

where q are the probability amplitudes. Then, the information entropy of the eigenvector 
\E > with respect to the basis set {\i >} is defined as 

S(E) = -J2\c l \ 2 \n\c l \ 2 . (15) 

i 

The idea is that, just as in the classical theory a dissolution of integrability (with the KAM 
mechanism) simply means the onset of chaotic motion, in quantum systems a dissolution 
of quantum numbers may indicate the onset of quantum chaos (see also [16]). In figure 
7 we show the information entropy S(E) of the eigenvectors of the Hamiltonian matrix of 
symmetry class (E,S) with respect to the axial symmetric basis set \J,k,S >, calculated 
for different values of the deformation parameter a. As expected, if the system has axial 
symmetry (a = \/2) then the information entropy S(E) is everywhere zero. By deforming 
the system, i.e. breaking the axial symmetry, S(E) becomes positive and it is larger in 
the central part of the energy spectrum. It is important to stress that in our system the 
dissolution of quantum numbers shown in Figure 7 does not have a classical analog because 
the classical Hamiltonian is always integrable. Thus, in our case, a large S(E) simply means 
a fully broken axial symmetry of the rigid rotor. 

CONCLUSIONS 

The main conclusion of this paper is that the asymmetric rotor is, like the harmonic 
oscillator, another pathological case with respect to the classical-quantum correspondence 
between integrability and Poisson statistics. In our opinion, the pathology of the asymmetric 
rotor model is more interesting because, unlike the harmonic oscillator, the asymmetric rotor 
satisfies the conditions of the Berry- Tabor theory. The presence of hidden symmetries could 
explain the pathological behavior of spectral statistics but such symmetries have not yet 
been identified. For the sake of completeness we remember that, as stressed by Rau [17], in 
classical mechanics, the asymmetric rotor and the nonlinear pendulum are intimately linked 
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and form the basis for many studies of nonlinear dynamics. Finally, we have shown that the 
information entropy of eivenvectors with the respect to the axial symmetric basis set gives a 
clear signature of the breaking of axial symmetry of the rigid rotator, but the rigid rotator 
is always classically integrable. 

V.R.M. is greately indebted to M.V. Berry and S. Gram for enlightening discussions. 
L.S. acknowledges M. Robnik and J.M.G. Gomez for fruitful conversations. 
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TABLES 



(E,S) (E,A) (0,S) (0,A) 
J even f + 1 { { j 
J odd ^ 

Table 1. Number of states in each sub matrix of the asymmetrical top Hamiltonian for 
a fixed J. 
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FIGURE CAPTIONS 



Figure 1: Density of levels p(E) of the four classes of symmetry with J = 1000. Pa- 
rameters: a = 1, b = a/2, c = a/5 and h = 1. 

Figure 2: Nearest neighbor spacing distribution P(s) of the four classes of symmetry 
with J = 1000. The dashed line is the Poisson prediction P(s) = exp (— s). Parameters: 
a = 1, b = a/2, c = a/5 and h — 1. 

Figure 3: Spectral rigidity A 3 (L) of the four classes of symmetry with J = 1000. The 
dashed line is the Poisson prediction A 3 (L) = L/15. Parameters: a = 1, 6 = a/2, c = a/5 
and 7l = 1. 

Figure 4: Density of levels p(E), Nearest neighbor spacing distribution P(s) and spectral 
rigidity A 3 (L) of the (E, S) class of symmetry, with J = 2000 (top) and J = 4000 (bottom). 
Dashed lines are Poisson predictions. Parameters: a = 1, b = a/2, c = a/5 and h — 1. 

Figure 5: Density of levels p(-E) of the (E 1 , 5) class of symmetry, with J = 1000. 
Parameters: b = a/2, c = a/5 and ft = 1. Different values of the deformation parameter: (a) 
a = a/2, (b) a = a/L9, (c) a = a/L5, (d) a = 1, (e) a = a/05, (f) a = 0. 

Figure 6: Nearest neighbor spacing distribution P(s) of the (E, S) class of symmetry, 
with J = 1000. Parameters: b = a/2, c = a/5 and U—l. Different values of the deformation 
parameter: (a) a = a/2, (b) a = a/L9, (c) a = a/L5, (d) a = 1, (e) a = a/o1>, (f) a = 0. 

Figure 7: Information entropy S(E) of the (E 1 , 5) class of symmetry, with J = 1000. 
Parameters: b = a/2, c = a/5 and U — l. Different values of the deformation parameter: (a) 
a = a/2, (b) a = a/L9, (c) a = a/L5, (d) a = 1, (e) a = a/05, (f) a = 0. 
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